set more off
estsimp logit forcea civila civilb  insta_xconst3 instb_xconst3 relpow jointcont2 jointally2 jointdem  instxconst3_civilb instxconst3_civila peaceyrs

setx civila 1 civilb 0  insta_xconst3 1 instb_xconst3 0 instxconst3_civilb 0 instxconst3_civila 1  relpow mean jointcont2 1 jointally2 0 jointdem 0 peaceyrs mean, 
simqi, prval(1) genpr(civilwarinstabhigh)
setx civila 0 civilb 0  insta_xconst3 1 instb_xconst3 0 instxconst3_civilb 0 instxconst3_civila 0 relpow mean jointcont2 1 jointally2 0 jointdem 0 peaceyrs mean
simqi, prval(1) genpr(instabhigh)
setx civila 1 civilb 0  insta_xconst3 0 instb_xconst3 0 instxconst3_civilb 0 instxconst3_civila 0  relpow mean jointcont2 1 jointally2 0 jointdem 0 peaceyrs mean
simqi, prval(1) genpr(civilwarhigh)
setx civila 0 civilb 0  insta_xconst3 0 instb_xconst3 0 instxconst3_civilb 0 instxconst3_civila 0 relpow mean jointcont2 1 jointally2 0 jointdem 0 peaceyrs mean, 
simqi, prval(1) genpr(civilwarinstablow)

gen fd_1=civilwarinstabhigh-instabhigh
gen fd_2=civilwarhigh-civilwarinstablow

centile fd_1, centile(2.5 97.5)

centile fd_2, centile(2.5 97.5)

gen sec_fd=fd_1-fd_2
centile sec_fd, centile(2.5 97.5)

drop b*
